SEPARATION OF VARIABLES AND BACKLUND TRANSFORMATIONS 
FOR THE SYMMETRIC LAGRANGE TOP 



VADIM B. KUZNETSOV, MATTEO PETRERA, AND ORLANDO RAGNISCO 

Abstract. We construct the 1- and 2-point integrable maps (Backlund transformations) 
for the symmetric Lagrange top. We show that the Lagrange top has the same algebraic 
Poisson structure that belongs to the sl{2) Gaudin magnet. The 2-point map leads to a real 
time-discretization of the continuous flow. Therefore, it provides an integrable numerical 
scheme for integrating the physical flow. We illustrate the construction by few pictures of 
the discrete flow calculated in MATLAB. 
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1. Introduction: the symmetric Lagrange top 

The Lagrange top is an integrable case of rotation of a rigid body around a fixed point in 
a homogeneous gravitational field, characterized by the following conditions [2, 4]: the rigid 
body is rotationally symmetric, i.e. two of its three principal moments of inertia coincide, 
and the fixed point lies on the axis of the rotational symmetry. A standard form of the 
corresponding equations of motion is given by the Euler-Poisson equations: 



;i.o.i) 



J = P X X, 
X = J X X. 
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Here J = {Ji, J2, Js) ^ M'^ is the vector of angular momentum of the body, P = (0, 0, a) G 
is the constant vector along the gravity field and x = {xi,X2, x^) e MP is the vector pointing 
from the fixed point to the center of mass. 

The symmetric Lagrange top is an integrable system with 2 degrees of freedom and the 
Hamiltonian 

(1.0.2) n = ^{J^ + Jl + J|) + axa, 

where Jk,Xk, k = 1,2,3 are six generators of the Lie-Poisson e(3) algebra defined by the 
following Poisson brackets: 

(1-0.3) {Jk,Jl}^Jm, {Jk,Xl}^Xm, {Xk,Xi}^0, 

{klm) is a cyclic permutation of (123). 

We will also use the complex conjugated variables J± — Ji ±iJ2 and x± — xi ±ix2, which 
have the brackets 

(1.0.4) {J3, J±} ^ TiJ±, { J+, J-} = -2iJ3, {J3,x±} ^ {xs,J±} ^^ix±, 
{ J+, = {x+, J_} = -2?:.T3, { J3, xs} = { J+, x+} = { J_, x^} = 0, 

{xk,xi} = 0, k,l = ±,3. 
The Lie-Poisson bracket (1.0.3) have two Casimir functions 

3 3 
(1.0.5) C, = J24, C2 = ^XkJk. 

k=l k=l 

Fixing their values one gets a generic symplectic leaf 

(1.0.6) Oei.c^ = {x, J\ Ci = Ci, C2 = C2}, 

which is a four-dimensional symplectic manifold. Hereafter we take Ci — 1 and C2 = i, 
corresponding to a unit vector x and a fixed projection £ of the angular momentum J on the 

vector X. 

Two commuting integrals of motion (Hamiltonians) of the symmetric Lagrange top are 
respectively, 

(1.0.7) n=-{J^ + Jl + Jl)+ax3 and J3, {n,J3} = 0. 

Obviously, the conservation of J3 is a direct consequence of the invariance under rotation 
about the direction of the gravity field. Using complex generators of e(3) we can write the 
Hamiltonian as 

(1.0.8) H = -{J+J- + J|) + ax3. 

2. From the sl{2) Gaudin magnet to the symmetric Lagrange top 

The symmetric Lagrange top can be derived from the s/(2) Gaudin magnet [5, 7], which has 
the 2x2 Lax matrix 

M.)^|:^(|_'y.«(j °)K;^11-1iri)' 

(2.0.2) Ag{u)^a + J2^, Bg{u) = J2^^, Cg{u) = Y,^ , 

j=l ■' J=l ■' 3=^ ■' 
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where Oj G C and a G M are the parameters of the model and m G C is the spectral 
parameter. The parameter a has the meaning of magnetic field's intensity. Local variables 
s|, s^, j = 1, are the generators of the direct sum of n sl{2) spins with the following 
Poisson brackets: 

(2.0.3) {s^ 4} = TiSjkst } = -'^iSjksl 

We denote the Casimir operators (spins) as Sj: 
(2.0.4) s^ = {s^r + spj. 

Fixing s/s we go to a symplectic leaf where the Poisson bracket is non-degenerate, so that 
the symplectic manifold is a collection of n spheres. 

The Lax matrix (2.0.1) satisfies the linear r- matrix Poisson algebra: 
(2.0.5) 

1 2 1 2 1 /l0\2/l0\ 

{Lg {u),Lg (v)} ^[r{u-v),Lg {u)+ Lg {v)l Lg^ Lg ^ ^1, Lg^ [q ^j^Lg, 
with the permutation matrix as the r matrix: 

/ 1 \ 

\ 1 y 

Equation (2.0.5) is equivalent to the following Poisson brackets for the rational functions 

Ag{u), Bg{u), Cg{u) (2.0.2): 

(2.0.7) {Ag{u),Ag{v)} = {Bg{u),Bg{v)} = {Cgiu),Cg{v)}=0, 

(2.0.8) {Ag{u),Bg{v)} = -^[Bg{v) - Bg{u)], 

(2.0.9) {Ag{u),Cg{v)} = -l-[Cg{u) - Cg{v)], 

(2.0.10) {Cg{u),Bg{v)} = -^[Ag{u) - Ag{v)]. 

The non-linear dynamics defined by the equations (1.0.1) is linearised on the Jacobian of 
the hyperelliptic spectral curve Fg of genus n — 1, 

(2.0.11) Fg : det{Lg{u) - v) = 0, 

which can be brought into the form 



(2.0.12) = AUu) + Bg{u)Cg{u) ^a'' + y( 

^-^ \U — tti 

The Hamiltonians Hj above are given by 



^1 



{u - ajY 



(2.0.13) H, = J2 '^'^''^'^''^'^''^ + 2asl 

These are integrals of motion of the s/(2) Gaudin magnet, which are Poisson commuting: 
(2.0.14) {Hj,Hk}^0 j,k^l,...,n. 
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Take a 2-site Gaudin spin chain for which n = 2 and the phase space is the direct sum of 
two sl{2) spins. We use the notations s] and s'j defined by sf = s^± is'j. Introduce a new 
matrix, called L{u), as follows: 

(2.0.15) L{u) = iL^^='\u)^^—L, + ^—L2 + a(' °. V 

where 

(2.0.16) L,^( /I , ^^+jf^)es«(2), k = l,2. 

Let us consider the following Lie algebras' isomorphism: 
(2.0.17) su(2)®su(2)^o(A). 

The Inonii-Wigner contraction from the rotation group 0{m + 1) to the Euclidean group 
E{m) [9], for m = 3, allows us to obtain the Lie-Poisson algebra e(3) for the classical 
Lagrange top, with six generators Jfe, x^, = ±, 3. 

Let us introduce the contraction parameter e G (0, 1]. Impose in the equation (2.0.15) the 
coalescence 02 = ai + e, e — > and redefine u — a\ = u. Then the Lax matrix reads 

(2.0.18) L(m) = i(Li + L2) + ^eL2 + a(^ Q ^.^+0{e^). 

In order to have a 2 x 2 Lax matrix for the Lagrange top we have to control that Li + L2 
and eL2 play the role of the following matrices respectively: 

or, in other words, we have to check the following Poisson morphisms: 
(2.0.20) sl + sl^J^, s+ + 4^J+, Sj: + S2"^J-, 

3 ~|~ — 

With a direct calculation of Poisson brackets it is possible to prove that morphisms (2.0.20) 
actually hold. 

Therefore we obtain the 2x2 Lax matrix for the Lagrange top 

We want to remark that in the Lagrange case the real parameter a describes the gravitational 
field's intensity. 

The spectral curve of the Lagrange top is the elliptic curve F : det(L(u) — v) — 0, 

2.0.22 r : -v^ ^ + + + ^ + a\ 

u-^ u 

The Lie algebra contraction o(4) — e(3) does not change the r-matrix structure of the 
model, that remains in the form (2.0.5)-(2.0.6) with the Lax matrix given by (2.0.21). 
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3. Separation of variables 

In this Section we construct the simplest separation of variables for the symmetric Lagrange 
top with the Lax matrix (2.0.21). The details of the approach can be found in [15, 10, 11]. 

The basic separation has only one pair {ui,Vi) of separation variables belonging to the 
spectral curve F (2.0.22). It corresponds to the standard normahsation vector ao = (1,0) 
and it is defined by the equations 

(3.0.1) (l,0)(L(^^i)-^;i)^ = 0, 

where (•)^ denotes the adjoint matrix. Equation (3.0.1) is the equation for the pole {ui,vi) 
of the Baker-Akhiezer function ^, which is defined as a properly normalised eigenfunction 
of the Lax matrix: 

(3.0.2) L(u)* = 'u*, {u,v)er, 

(3.0.3) ao • * = 1. 

It is easy to see that the equation (3.0.1) gives the following separation variables: 
(3.0.4) B{ui) = 0, vi^ -A{ui). 

Explicitly, the pair of the canonical separation variables is 

(3.0.5) ui^-^, vi^i—(j3-—X3 

In order to define the map between the initial and separation variables, one has to add 
an extra pair of canonical variables. This pair of variables is taken from the asymptotics 
{u — > oo) of the elements A{u) and B{u) of the Lax matrix. We obtain 

(3.0.6) J-, -i^. 

Indeed, using the Lie-Poisson algebra e(3) (1.0.4), it is easy to check that the new variables 
are canonical: 

(3.0.7) {uk,ui} = {vk,vi} ^ 0, {vk,ui}^5ku k,l^ 1,2. 

Remark that the new separation variables (mi, «2, ^^i, ^2) are defined in a complex domain. 
We can write down the complex generators of e(3) in terms of the separation variables: 

(3.0.8) J_ = U2, 

1 ulvf 



J+ = H V2U2 2" 

U1U2 UIU2 U2 

J3 = iV2U2, 
X_ = -U1U2, 

1 2 ^2 f4'<^l 
X+ — V2U1U2 + 2V1V2U1 . 

U1U2 U2 

= iuiiuiVi - U2V2) . 



The separation equations are: 

/ (. 1 

(3.0.9) T-i ^ u\\ iavi - ] - 7-^ - aJsUi, 

\ 2 J ui 2u\ 

(3.0.10) J3 = iv2U2. 
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One can use the above separation of variables to integrate the model in terms of elliptic 
functions. Although, because the separation variables are complex, it is not easy to apply the 
reality conditions. In the next Section we describe the standard integration of the symmetric 
Lagrange top. 

We finish this Section by presenting the canonical separating transform through a gener- 
ating function. 

Let us fix a special representation of the e(3) algebra in terms of the Darboux coordinates 
Qj^PjJ = 1,2: 

(3.0.11) J_ = qi, 

J+ = qipl + 2q2PiP2 - 2i£pi - 2ip2, 

Js = -i{qiPi + ?2P2) - ^, 

X- = q2, 

x+ = q2pl-2ipi, 

X3 = -iq2Pi - 1, 

(3.0.12) {q^^qi} = {p^^pi} = 0^ {qk,Pi} = Sm, A:,/ = 1,2. 

Notice that in this representation the variables J_ and a;_ do not depend on the momenta 
and the variables J3 and xs are linear in the momenta. 

Finally, the canonical transformation {qi,q2,Pi,P2) — ^ (wi, ^2, t^i, 1^2) is defined by the 
generating function 

(3.0.13) F{qi,q2 | Vi,V2) = i — + ii\nq2 - qiV2 + , 

q2 qi 

that is 

/Q n dF{qi,q2 \ Vi,V2) dF{qi,q2 \ Vi,V2) 

(3.0.14) = — , Ui^ — , 1 = 1,2. 

4. Integrating the model 

The standard integration of the Lagrange top is performed through Eulerian variables which 
allow the reduction to one degree of freedom by using the first order integral J3 to separate 
out the angle tp and its conjugated momentum p^ [1]. 

Let us introduce the Eulerian angles {6,(p), 6 & [Ojh], if E [0,27r), to parametrize the 
unit sphere xf + xl + = 1, which together with the conjugated momenta {p0,Pyp) give a 
representation of the e(3) generators in terms of Darboux coordinates: 

/ . r, -, \ T £-p^cos9 . . . _ 

(4.0.1) Ji = po cos (p-\ sm</7, Xi = sm (/? sm t/, 

sin 6' 

J . e-p^cos9 

J2 = —pe sm ip H cos ip, X2 = cos (p sm 6, 

J3 = p^, X3 = cos9, 

(4.0.2) {pe,e} = {p^,ip} = l. 

One immediately arrives at two one- dimensional equations 

pi-2ip^cose + e 

(4.0.3) H = ^ + ^^^ f^^- + acos^, 

2 2 sm U 

(4.0.4) J3 = p^. 
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which lead to integration of the model in terms of the Weierstrass elliptic function p{x). The 
time evolution of the angles and of their conjugated momenta is the following: 



In this paper, following the approach of [13, 14], we look at the Backlund transformations 
(BT's) for finite-dimensional (Liouville) integrable systems as special canonical transforma- 
tions, thereby taking a Hamiltonian point of view. Such BT's are defined as symplectic, 
or more generally Poisson, integrable maps which are explicit maps (rather than implicit 
multivalued correspondences) and which can be viewed as time discretizations of particular 
continuous flows. 

The most characteristic properties of such maps are: (i) a BT preserves the same set of 
integrals of motion as does the continuous fiow which it discretizes, (ii) it depends on a 
Backlund parameter A that specifies the corresponding shift on a Jacobian or on a gener- 
alized Jacobian [14, 3], and (iii) a spectrality property holds with respect to A and to the 
conjugate variable /i, which means that the point (A, /S) belongs to the spectral curve [13, 14]. 
Explicitness makes these maps purely iterative, while the importance of the parameter A is 
that it allows for an adjustable discrete time step. The spectrahty property is related with 
the simplecticity of the map [14]. 

In this paper wc construct BT's for the symmetric Lagrange top starting from the results 
obtained in [8]. As shown in previous sections, the Lagrange top has the same algebraic 
Poisson structure that belongs to the Gaudin magnet. This allows to choose the same 
ansatz for the matrix M{u) that has been used in [8]. 

We have here to remark that an elegant alternative approach to integrable time discretiza- 
tions of continuous hamiltonian fiows has been carried out by Yu. B. Suris and A.I. Bobenko 
[6, 16]. In particular, in [6] they construct a discrete time Lagrange top: notice howewer that 
their discrete Lax matrix is deformed with respect to the continuous one, which is recovered 
in the limit when the time step goes to zero. 

5.1. One-point BT. A (one-point) Backlund transformation for the Lagrange top is equiv- 
alent to the following similarity transform on the Lax matrix L[u): 



with some non-degenerate 2x2 matrix M{u, A), simply because a BT should preserve the 
spectrum of L{u). The parameter A G C is called a Backlund parameter of the transforma- 
tion. Let us introduce new ~-notations for the updated variables: 



(4.0.5) 




5. Backlund transformations (BT's) 



(5.1.1) 



L{u) I — > M{u; X)L{u)M-\u, A) Vm, 



A e C, 



(5.1.2) 
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We are looking for a Poisson map that intertwines two Lax matrices L{u) and L{u): 

(5.1.3) M(u;X)L{u) ^L(u)M(u;X) ^u. 
Let us take 

(5.1.4) M(«;A) = ^ ^ , detM{u;X) = u-X. 

Let us stress that the number of zeros of det M is the number of essential Backlund param- 
eters. Here the variables p and q are indeterminate dynamical variables. 

The ansatz (5.1.4) for the matrix M came from the simplest L-operator of the quadratic 
r-matrix algebra 

(5.1.5) {L iu),L (v)} = [r{u - v), L (u) L (v)], 

with the same r-matrix (2.0.6). Notice that we have a peculiar algebraic situation, as in 
Gaudin models: a Lax matrix which satisfies a linear r-matrix algebra requires a M matrix 
which comes from a quadratic r-matrix algebra. Recall that in the Toda lattice [13] and in 
the DST model [12] both L and M are derived from the same quadratic r-matrix algebra. 
This fact shall show up in the construction of the quantum analogue of these Backlund 
transformations . 

Comparing the asymptotics in — > oo in both sides of (5.1.3) we readily get 

(5.1.6) Ja^Ja, g = ^ . 

If we want an explicit single- valued map from L{u) to L{u) we must express M{u, X) , and 
therefore p and q, in term of the old variables. To solve this problem we use the spectrality 
of the BT. As well as the equation (5.1.3) that our map satisfies, it will be parametrized by 
a point P = (A, /x) G T. Notice that there are two points on F, P = (A, and Q = (A, — //), 
corresponding to the same A and sitting one above the other because of the elliptic involution: 

(5.1.7) (A,/x)Gr: det(L(A) - /i) = ^ + det{L{X)) = 0. 

As shown in [8], this spectrality property, used as a new datum, produces the formula 

^" B(X) - A(A)+/.' 
where A and // are bound by the equation for the elliptic curve 

2 f 1 2e 2n 2J3 2 



(6.1.9) = + _ + _ + + a 

Now the equation (5.1.3) gives an integrable Poisson map from L{u) to L{u). The map is 
parametrized by one point (A, //) G F. Explicitly, it reads 

(5.1.10) Ja = Ja, 

J_ = X- + J-{pq - X) - 2pJ3, 

J+ = 2qa, 

X3 = X3 + pJ+ - gx_ - qJ-{pq - X) + 2pqJ3, 

X- = x_ {2pq - A) - 2pxz -p'^J+ + J-pq{pq - A) - 2p^g Ja, 

x+ = J+ - J^-{pq - X) + 2qJs. 
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If we refer to the real generators of e(3) the map takes the following explicit form, 



(5.1.11) 



Jl 



J2 



2 2 4q; 

i 

+7r{J2J3 - qJiJ2 + iJiJ^), 

la 

-iaq + Uxi + ix2) - ^(Ji + iJ2) + t^{Ji - Jl) + 

+—(72/3 - gJi J2 - iJiJz), 
la 

J3, 



Xi 



1 X 
(^Q^ + -jjyJi + iJ^) - ^(-^1 - iJ2) - 2 (^1 ~ ^^2) + qJi + 



2a 
8a^ 



A 



2^3 



-{Jl - J{) + (Jl -iJ2)[xi - — ] - ix2{Ji + iJ2) + ^AJl J2 



+ 



X2 



qJl{Ji - 3iJ2) + qJl{J2 - 3iJi) - '(Ji ' J2)\Ji + J2) - 2Js{Ji - iJz) 
iq'^ 



-iaqX + ^(Ji - zJa) - ^(<^i + ^-^2) - y (a^i - ^a;2) - iqJs + 



iq 
^2a 



A 



2^3 



-( J2 - Jl) + (Jl - iJ2) \ xx ia;2( Jl + iJ2) + iAJi J2 



+ 



iq 



8a2 



gj2( Jl - 3iJ2) + gJ2'( J2 - 3iJi) - -( Jl - J2)'( Jl + J2) - 2 J3( Jl - iJ2f 



1 

+- 
a 



^-{Jl + Jl) - %-{Jl - Jl) - iq{J2J3 - gJi J2 + ^Jl J3) 



2^3 = 2:3 - q{xi - 1x2) + gA( Jl - 2 J2) + 

"1 -2 

■ 2 

where we have used (5.1.6) to express p in terms of the old variables. 

Notice that the above one-point BT is a complex map, so it is a non-physical Backlund 
transformation . 

5.2. Symplecticity of the one-point BT. Wc give a simple proof of symplecticity of the 
constructed map by finding an explicit generating function of the corresponding canonical 
transformation from the old to new variables. 

First, because the Casimir functions do not change under the map. 



(5.2.1) 
(5.2.2) 



ly I nf rp rp^ I rp rp 1 

^(J+x_ -\- J-X+) + J3X3 = ^{J+x- + J-X+) + J3X3 = £, 



we can exclude the variables x+, J+ and J_, using the following substitutions: 



(5.2.3) 
(5.2.4) 



1-4 

X- 



x+ 



J 

J+ — — 

X- 



2J3a;3 J 



2J3J3 



2 (1 2:3) , 



x+ 



^(1-sl). 
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Now we have only four (old and new) independent variables: x^,X3, J_, J3 and x+, X3, J+, J3. 
We write the 1-point BT as a canonical transformation defined by the generating function 
Fx(a;_, J_ I 5;+, J+) written in terms of x_, J_ and 5;+, J+: 

dFx{x-, J-\x+, J+) 



(5.2.5) X3 — ix- 



dJ. 



. (9Fa(x_, J_ I x+, J+) aFA(a;_, J_ I J+) 

J3 = IX- h iJ- 



X^ — IXj^ 



dx- dJ- 
dFx{x_ , J_ I x+, J+) 



~ _ dFx{x_,J_ \ x+,J+) ~ dFx{x-,J-\x+,J^ 

Jr — IX+ — h iJj, 



dx+ dJ+ 
With the help of (5.1.6) we rewrite equations (5.1.10) of the map in the form 

(5.2.6) X3 = + k, fc^ = 1 + Xx_x+, 

2a 

It is now easy to check that the function 

(5.2.7) Fx{x_,J_ I J+) = -^ - ^J + ^ - i ) + z£log - zaA, 

2q; \x-X-^Xj 1 — k 

solves the equations (5.2.5)-(5.2.6). Symplecticity of the map is therefore proven. 

Alternatively, we can derive (5.2.7) directly from the generating function constructed for 
the sl{2) Gaudin magnet in [8]. The one-point BT in that case is given by 

(5.2.8) - is^ — , - isj g— , 
where 

(5.2.9) Fx{s^,...,s- I S'i+,...,s„+) = --^ V s/Sfc^-^y" ( 2zj + 5^ log - ia\ 

(5.2.10) = s]-{a,-\)s,+sj, s] = {s]f + s- = (S^f + sfs+, j = l,...,n. 

Set n = 2 in (5.2.9) and impose the coalescence a2 = Oi + e, e — > 0, similar to the derivation 
of the Lax matrix for the Lagrange top in Section 2. Now redefine A — ai = A and obtain 
(5.2.11) 

^A(5r,52 I 5i+,S2+) --^Y. ^7^-^^ f20 + s,log^^^') - I'^e - laX + 0{^) , 

where 

(5.2.12) (^^ = s^. + Xsj+s-, j = l,2. 
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(5.2.13) Cl = sl + X{j-J+ 



t2 T T 2£ 1 
J| + J+ J_ ^ — 



(5.2.14) 



Ci = sl{l + Xx-x+), 



Now it is easy to check that (5.2.11) with (5.2.13)-(5.2.14) coincides with the generating 
function (5.2.7). 

5.3. Spectrality of the one-point BT. The spectrality property of a Backlund transfor- 
mation [13] means that the two components, A and /i, of the point (A, //) G F parametrizing 
the map are conjugated variables, in the sense that 

dF 



(5.3.1) 



where F is the generating function of the BT. 

We now show the spectrality property for the one-point BT constructed in the previous 
section. Using the equations (5.1.6), (5.1.8), (5.2.6) and (5.2.7) we obtain 



(5.3.2) fx 



1 

A2 



A 



J_ J+ 1 , 
— + \ + ak 

a:_ x+ A 



^(A)-^i.(A)4 

g^A(.r-.-/-|7+.J+) 
dX 

5.4. One-point BT in Eulerian variables. As shown in Section 4 one can represent the 
Lie-Poisson algebra e(3) using the canonical realization (4.0.1). Now we rewrite the one-point 
BT canonical map in terms of these variables 

(5.4.1) Bx: (6l,v?,pe,p^) I — > {d,lp,pe,p^), 

where tilde refers to the updated versions of the generators 



(5.4.2) 



~ ~ ~ i — Pu> COS 9 _ 
Ji = pe COS (fi H — sm (p, 

sin 9 

~ ~ i — Pw cos 9 _ 
Pe sin ip H — cos (p, 



J2 = 
J3 = P<p 



sin^ 



Xi — sin (/7sin 9, 

X2 — cos (f sin 9, 

X3 = cos 9, 



(5.4.3) 



{pe,9} = {p<p,(p} = 1. 
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Recall that we have fixed the values of the two Casimir functions as follows: Yl 
and Ylfk=i-^>'Jk — ^- Using the equations (5.2.6) we can write down the one-point BT as: 



k=l 



(5.4.4) 



(5.4.5) 



with 
(5.4.6) 



= 

P<P = 
Pe = 

P9 = 



1 
k 

Pip, 



\a{cos6 + cos 6') + 



P-1 
2A 



+ 2a\ 



(cos^- l)(cos^- 1) 



sin^ 

i 

sm.9 



{£ - cos 6) — 
{p^ cos 9 - 



+ 



2iaX sin ^(cos 9 — k) 

2ia\ sin ^(cos 9 — k) 
k'^-1 



k'^ — X(f)(f) sin 9 sin 9, 



4> 



Remark that in the last formula we have introduced two new variables (j) ^iid </'- it is a 
convenient choice since the map depends only on e*'^ and e^*'''. 

We want to focus our attention upon the symmetry of the canonical transform (5.4.4)- 
(5.4.5), which is symmetric under the exchange {9^Lp) ^ {—9,—(p). Notice that the 0,0 
dependence is trivial: it involves only the product 00 (i.e. (f — lp), which is symmetric under 
the change <-> 0. The dependence upon 9, 9 is more interesting. Let us introduce the 
exchange operator Vg g such that 

(5.4.7) V,~,f{9,9)^f{l9l Vl^^l, 

where f{9,9) is a generic function. For instance we notice that Vggk = k, where k is the 
function in (5.4.6). If we consider formulae (5.4.4)-(5.4.5) it is easy to notice that 

(5.4.8) 'PggP^=P^, VggP0 = -p0. 

Our aim is to write down the canonical transform in the following form: 

dF^{9,cl)\0,^) ^ dFx{9,cj)\9,^) 



(5.4.9) 



Pe 



Pv = 



89 



Pe 



89 



84> 



^ ~ 8Fxi9,<p\9,4>) 

P^ = Z0 - 



where Fx{9 , 4)\9 , (p) is the generating function of the one-point BT. Since equations (5.4.8) 
hold the immediate consequence is that 

(5.4.10) F,{9, (P\9, 0) = Vg gF,{9, 0|^, 0) = Fx{9, ^\9, 0). 
It is possible to prove that the generating function is: 

(5.4.11) Fx{9, 4>\9, 0) = -taX{2 In 00 + 1) + Ax{9, 4>\9, 0) + r,~,Ax{9, (f>\9, 0), 
where 

(5.4.12) Ax{9, 01^, 0) = ao + a,E{^\rj) + a2F{C\v) + asU{C; 
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with 



ao 



(^Utt) -^ix + ¥^^^ (S) + 2^"^ |ln(sin^)(l + cos^) - k+ 

+F5f cos^ + ^(1 + cos ^ cos ^) + ^J:^:^,^ + In (^) + 
+ |cos^ln(l±|) + 3im±i)}, 



(5.4.13) < 



ai 



0.1 — 



as 



2iaS 
kcosB 



2iaS 
kcosO 



:i + cos^) (3A-^^,)- 



A</.202(l-cos6i) 



fe cos a sm p 



(5.4.14) 



r m\v) 

I n(C;^|r7) 



(5.4.15) 



Jq^(1 — sin^ ?7sin^ z)^''^d2;, 

(1 - sin^ ?7sm^ ^)"^/^dz, 
/o^(1 - C^sin^ '2)~^(1 - sin^rysin^ z)~'^/'^dz, 



k 



I C 



sin0 
fc2-l+sine 



fc2-l+sin6> 
■fe2-l-sin6» ' 



sin 6 



In the zero-field case {a — 0) the generating function simply reads: 



(5.4.16) 



.(a=0). 



a In 



1 + k 



sin ^ sin 9 



k 
A 



l-^(cose + l)(cose + l) 

5.5. Two-point BT. According to [8], we now construct a composite map which is a prod- 
uct of the map Bp^ = B(^Xi,ixi) and Bq^ = B^X2-H2)' 

(5.5.1) Bp,,Q, = Bq., o Bp, : L{u) ^ L{u) ^1 (u). 

Two maps are inverse to each other when Ai = A2 and jii = Hi- This two-point BT for the 
Lagrange top is defined by the following 'discrete-time' Lax equation: 

(5.5.2) M{u]\i,\2)L{u) =L{u)M{u]\u\2) V-u, Ai,A2eC, 
where the matrix M(m; Ai, A2) is 

u — \i + xX X 

—x^X + (Ai — A2)x u — \2 — xX 



(5.5.3) 
(5.5.4) 



M{u; Ai,A2) 

det M{u; Ai , A2) = (w - Ai) {u - A2) 
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Spectrality property with respect to two fixed points (Ai,//i) e F and {X2, fJ'2) G T give 
A(Ai)-//i C(Ai 3(A2)-//2 5(A2) J+ 



(5.5.5) x = 



(A2-Ai)S(Ai)5(A2) 



^ ■ ■ ^ 5(Ai)(A(A2) + ^2) - 5(A2)(A(Ai) - ^1) 

(Ai-A2)(A(Ai) + /xi)(A(A2)-/x2) 
(A(Ai) + ^i)C(A2) - (A(A2) - ^2)C(Ai) 

(A2 - Ai) B (Ai) 5 (A2) 

B (A2)(l (Ai) + Ail)- S (Ai)(3 (A2) - fi2) 

(Ai - A2)(l (Ai) - I^i)(a (A2) + At2) 

(1 (A2) + At2) 5 (Ai) - (1 (Ai) - III) C (A2) 

J-- J- 
2a 

Now we have two Backlund parameters Ai, A2 G C. The above formulae give several equiv- 
alent expressions for the variables x and X since the points (Ai,/^i) and (A2,//2) belong to 
the spectral curve F, i.e., are bound by the following relations 

(5.5.7) 111 = A\\k) + B{\k)C{\k) (Afc)+ B (A,) 5 (A,), = 1, 2. 

Together with (5.5.5)-(5.5.6), the formula (5.5.2) give an exphcit two-point Poisson inte- 

grable map from L{u) to L {u) (as well as its inverse, i.e the map from Z {u) to L{u)). The 
map is parametrized by two points Bp^ = and Bq^ = B^x^^^^^y Let us introduce the 

notation 

(5.5.8) AA = Ai - A2, Ao = . 
The map exphcitly reads: 

(5.5.9) J3=^3, J-=J--2X, J+= J+ + 2x{AX-xX), 
xs = X3- a(AX- xX)J^ + XJ+ + 2xX(AX- xX), 

x_ = x- + {2xX + AX)J--2XJ3-2x(^xX + Xo-^ 

{2xX + AA) J+ + 2a;(AA - xX) J3 - 2x(AA - xX) (^xX - Aq - . 
It is useful to write down the Backlund map in terms of the real generators of e(3): 

(5.5.10) Ji^ Ji + xAX-X(l + x^), J2^ J2-i[xAX + X{l-x'')], Js^ J3, 
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xi = xi-iJ2(2xX - AX) + JslxAX- X(l + x'^)] - 



-\xX\l + x^) + X 



An- 



AA 



-xUXo + 



3AX \ 



xAX An + 



X2 — 



X2 + iJi{2xX - AX) - iJs [xAX + X{1 - x^)] 
-i\xX^{l-x'^)+X 



AX 2 A 3AAA 



AA 



xAX ( Ao + 



X3 = X3-Ji[xAX-X{l + x^)]+iJ2[xAX + X{l-x^)]-2xX{xX-AX). 

Obviously, when Ai = A2 (and fj,i — the map turns into an identity map. With a direct 
calculation one can show that the map (5.5.10) sends real variables to real variables provided 

(5.5.11) Ai = A2 = A e C. 

Therefore, the two-point map leads to a physical Backlund transformation with two real 
parameters. 

5.6. Two-point BT as a discrete-time map. The two point BT constructed above is 
a one-parameter (Ai) time discretization of a family of flows parametrized by the point 
Q2 = (A2,— At2), with the difference AA playing the role of a time step. Recall that the 
physical time step is iAA since AA e iR. 
Consider the limit 

(5.6.1) zAA = i(Ai - A2) = e, e ^ 0. 

It is easy to see from the formulae of the previous section that 

AiX2) /i2 C(A2) 



(5.6.2) 



X = Xo + O{e), 



Xq 



X 



eXo + 0(e2), Xo = 



The matrix M has the following asymptotics: 



(5.6.3) M{u;X2)^iu-X2) 1 + - , 

2n2{u - X2 

If we define the time derivative L{u) as 

L (u) - L{u) 



B{X2) A(A2)+//2' 

iB{X2) 

2/12 

A{X2)+fi2 B{X2) 

C{X2) -A{X2)+I12 



+ 0{e'). 



(5.6.4) 



L{u) = lim ■ 



then in this limit we obtain from the equation M{u)L{u) —L {u)M{u) the Lax equation for 
the corresponding continuous flow that our BT discretizes, namely: 

i:(A2) 



(5.6.5) 



L{u) 



2//2 



(u - A2 



-,L{u) 



This is a Hamiltonian flow with ifj,2, 
\/=A(A^ 



(5.6.6) 



IH2 



\^/% \ /I 2i 27Y 2J3 „ 

i?(A2)C(A2) = WT4+T3+Tr + T^ + «'' 

Ao An A9 A2 



as the Hamiltonian function, 
(5.6.7) 



Lij{u) = i{fj,2,Lij}. 
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Hence, the constructed two-point BT discretizes a one-parameter family of flows, labelled 
by the arbitrary value of A2. Now, let us consider the following limit: 

(5.6.8) A2 = iv, V ^ 0. 
In this limit we obtain the following expression for ^2' 

(5.6.9) = + 
The Lax equation (5.6.5) turns into: 



(5.6.10) L{u) = ^ 



1 , L(u) 



J-X^—J^X- J3X-—J-X3 I ax- 

J+xs—Jsx^ ax+ J-\.x-—J—X-\- 



„2 



2v? 

which is the flow with the Hamiltonian Ti (1.0.2): Lij{u) = {Ti, Lij{u)}. 

It is possible to show that one iteration of the constructed two-point BT is equivalent to 
moving from time to time 1 along the flow with the Hamiltonian (interpolating flow): 

(5.6.11) H{Xi, A2) = i / ' n{X)dX. 

5.7. Numerics. Below is the MATLAB program, which shows the real 2-point map from 
the Section 5.5. 

function bt3 = l(aa,lambdar,lambdai,N, JJl, JJ2, JJ3,xxl,xx2,xx3) 
alpha=double(aa) ; 

lainbdal=double (complex (double (lambdar) ,double(lambdai) )) ; 
lainbda2=double (conj (lambdal)) ; 
J3(l)=double(JJ3) ; 

Jp(l)=double(complex(double(JJl) , double (JJ2))) ; 
Jm ( 1) =double ( con j ( Jp ( 1 ) ) ) ; 
x3(l)=double(xx3) ; 

xp(l)=double(complex(double(xxl) , double (xx2))) ; 
xm(l)=double(conj (xp(l))) ; 
Cl=double(x3(l)*2+xp(l)*xm(l)) ; 

C2=double (xp ( 1) * Jm ( 1) +xm ( 1) * Jp (1 ) +2*x3 ( 1) * J3 ( 1) ) ; 
H =double(J3(l)-2+Jp(l)*Jm(l)+2*alpha*x3(i)) ; 

mul=double(sqrt (double (alpha*2+2*alpha*J3(l)/lambdal+H/lambdal*2+C2/lambdal"3+Cl/lambdal"4))) ; 
mu2=double (conj (mul) ) ; 
L0= [alpha ; -alpha] ; 

LJ(1,1,1)=J3(1); LJ(l,2,l)=Jm(l); LJ(2, 1 , l)=Jp(l) ; LJ(2,2, 1)=-J3(1) ; 
Lx(l,l,l)=x3(l); Lx(l,2,l)=xm(l) ; Lx(2, 1 , l)=xp(l) ; Lx(2,2, l)=-x3(l) ; 
for k=2:N 

x=double(((alpha-mul)*lambdal"2+J3(k-l)*lainbdal+x3(k-l))/(Jm(k-l)*lambdal+xm(k-l))) ; 

ddl=double((Jm(k-l)*lambdal+xm(k-l))*((alpha+mu2)*lambda2-2+J3(k-l)*laiiibda2+x3(k-l))) ; 
dd2=double((Jm(k-l)*lambda2+xm(k-l))*((alpha-mul)*lambdal-2+J3(k-l)*lainbdal+x3(k-l))) ; 

X=double (double ( (lambda2-lambdal) * (Jm(k-l) *lambdal+xm(k-l) )* (Jm(k-l) *lambda2+xm(k-l) )) /double (ddl-dd2) ) ; 
M0= [double (-laiiibdal+x*X) X; 

double (-x"2*X+(lambdal-lambda2)*x) double (-lambda2-x*X)] ; 
LJ(: , : ,k)=double(LJ(: , : ,k-l)+MO*LO-LO*MO) ; 
Lx(: , : ,k)=double(Lx(: , : ,k-l)+MO*LJ( : , : ,k-l)-LJ(: , : ,k)*MO) ; 
J3(k)=LJ(l,l,k); x3(k)=Lx(l , 1 ,k) ; Jm(k)=LJ(l,2,k) ; xm(k)=Lx(l,2,k) ; 

end 

bt3=plot3(real(xm) ,real(l*xm) ,real(x3) , ' :oy' , 'LineWidth' ,0.5, 'MarkerEdgeColor ' , 
' b ' , ' MarkerSize ',1.5,' MarkerFaceColor ' , ' b ' ) ; 
axis ('equal') 

The input parameters are: 

• aa = a, 

• lambdar = JJA, 

• lambdai= ^X, 

• N — number of iterations of the map. 
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• JJl,JJ2,JJ3,xxl,xx2,xx3 = initial values of Ji, J2, J3, xi, 0:2, 0:3. 
The output is a 3D plot of consequent (blue) points {xi, X2, X3) linked by the yellow lines. 
Notice that all these points lie on the sphere xl+xl + x"^ = constant, of some radius defined 
by the initial data. 




Figure 1. (a;9iA,OA;V;Ji,J2,J3,a;i,a;2,3;3)=(3.345;0,0.02;1000;2.34,6.4,8,1.219,-0.78,6.77) 



6. Concluding remarks 

We have constructed the 1- and 2-point Backlund transformations for the symmetric 
Lagrange top following the approach of [13]. The application of the constructed maps as 
exact numerical integrators of the continuous flows is considered. As shown in Section 2 
the Lagrange top has the same linear algebraic Poisson structure that belongs to the sl{2) 
Gaudin magnet. This allowed us to choose the same ansatz for the matrix M{u) (5.1.4) that 
has been used in [8]. 

Notice that we had the following algebraic situation, as in Gaudin models: a Lax matrix 
(5.1.2) which satisfies a /mear r- matrix algebra (2.0.5) requires a M matrix which comes from 
a quadratic r-matrix algebra (5.1.5). Recall that in the Toda lattice [13] and in the DST 
model [12] both L and M arc derived from the same quadratic r-matrix algebra (5.1.5). 
This fact shall show up in the construction of the quantum analogue of these Backlund 
transformations. Indeed an interesting problem is the quantization of the constructed maps, 
namely the problem of Baxter Q-operators for the symmetric Lagrange top. Actually this 
problem leads to the study of product formulae for Heun functions. This work is in progress 
and the results will be reported in a separate paper. 
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Figure 2 . (a;g?A,3!A;Ar; Ji, J2, J3,3;i,a;2,3;3)={l;500,5000;4000;3.4,-4.6,-6.2,-2.19,0.89,0.77) 
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